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FROM TRANSIENT DATA 
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SUMMARY 

Methods are presented that use general correlative time- 
response input and output data for a linear system to determine 
the frequency-response function of that system. These methods 
give an exact description of any linear system for which such 
transient data are available. 

Examples are shown of application of a method to both an 
underdamped and a critically damped exact second-order system, 
and to an exact first-order system with and without dead time. 
Experimental data for a turbine-propeller engine showing the 
response of engine speed to change in propeller-blade angle are 
presented and analyzed. 

INTRODUCTION 

A basic problem confronting the control designer is that 
of determining the behavior of the controlled system under 
varying conditions of operation. This question becomes 
particularly acute when various system parameters, for 
example, engine shaft torque or turbine-inlet temperature, 
must be closely controlled in order to prevent system damage 
or even failure. The problem then becomes one of matching 
transient behavior of the control to that of the uncontrolled 
system in order to attain a desired response of the controlled 
system. 

General methods do not exist for the solution of the 
equations of motion of nonlinear systems. As a result, the 
analysis of such systems may be impracticably difficult, if, 
indeed, a solution can be found at all. Because the behavior 
of many nonlinear systems may be satisfactorily approxi- 
mated by the assumption of system linearity and because 
general mathematical methods and techniques for the analy- 
sis of linear systems are readily available and (in comparison 
with present nonlinear methods) relatively simple to apply, 
the assumption is generally made that the system being 
studied is linear. The methods of this report are based on 
such an assumption. 

Approaches utilized in the analysis of linear systems 
(reference 1, pp. 17-18) are: (1) transient analysis, in which 
the characteristic time response of the system is determined 
for standard inputs such as step and impulse functions, or 
(2) frequency analysis, in which system behavior is char- 
acterized by the steady-state amplitude and phase relations 
of the system input and output for sinusoidal inputs of 
various frequencies. 

For inputs such as step or impulse functions, the inherent 
system characteristics might be obtained by fitting an 


equation to the output function and from it deriving the 
differential equation of the system. If the input and output 
functions were of any general form, fitting differential 
equations to the data might still be possible, but for systems 
of inherently high order the accuracy of such a procedure 
would be low. In general, use of the time-response form of 
system description involves dealing with convolution in- 
tegrals (reference 2, p. 54) when the unit so described is to 
be combined with other units. For complex systems, this 
descriptive form is not readily adapted to manipulation. 

A linear system is known to be characterized by its 
steady-state response to all frequencies of sinusoidal inputs. 
The frequency-response form is useful for the description of 
linear systems because of the ease with which various system 
characteristics can be manipulated in dealing with combina- 
tions of units. For example, the over-all amplitude ratio of 
output to input for any frequency of input- to a system con- 
sisting of several units in series is obtained by a simple 
multiplication of the amplitude ratios for the several units. 
A possible difficulty in the determination of the steady-state 
frequency response is the necessity for maintaining a sinu- 
soidal input of constant amplitude and frequency for a 
length of time sufficient to insure disappearance of transient 
effects in the output. 

The nature of real physical systems may make actual 
imposition or measurement of step or impulse inputs im- 
possible; in addition, sinusoidal inputs may prove impracti- 
cable. The need thus arises for feasible methods capable of 
handling data of any general form. Such methods that use 
general correlative time-response data for system input and 
output to deter min e the frequency-response characteristics 
of the system were developed at the XACA Lewis laboratory 
in 1948-49 and are presented in this report. These methods 
give an exact description of any linear system for which such 
data are available. 

Three exact methods of obtaining the frequency-response 
function of a system are shown for system-equilibrium final 
conditions. A graphical approximation to one of the 
methods is also given. Modification of one of the exact 
methods to account for oscillatory final conditions is shown 
and treatment of dead time is presented. As an illustration 
of one of the methods, three examples are given, two based 
on analytically determined system-response curves and one 
on experimental data obtained at this laboratory for a 
turbine-propeller engine. 
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SYMBOLS 

The following symbols are used in this report: 


A 

amplitude of system steady oscillation 

a,b 

constants 

F(io>) 

system frequency-response function 

F{p) 

system transfer function 

m 

general function of time 

G{p) 

system transfer function including dead time 

9(t) 

polynomial approximation to y{t) 

hit) 

normalized Gaussian error distribution or prob- 
ability pulse 

j 

order of derivative of u(t) 

L 

Laplace transform 

P 

complex number 

2 

degree plus one of g k (t) 

R 

amplitude ratio of F(i<a) 

t 

time, seconds 


unit step function 

u'(t) 

unit impulse function 

X 

general time-dependent input 

y 

general time-dependent output 

P 

frequency of system steady oscillation, radians 
per second 

At 

system dead time, seconds 

f 

damping ratio 

e 

phase angle of F(iu), radians 

a 

standard deviation of Gaussian curve 

T 

system time constant, seconds 


phase angle of system steady oscillation, radians 

U 

frequency of system sinusoidal input, radians 
per second 

Subscripts: 

0 

last term of differential equation 

/ 

final value 

k 

general term of summation 

m 

highest order derivative in expression for input 

n 

highest order derivative in expression for output 

r 

next to last term of summation 

Superscripts: 

m,« 

first, second, third, mth, and nth time deriva- 
tives, respectively 
ANALYSIS 

The system considered in the following derivation is 

assumed to be linear. In particular, behavior of the system 

is assumed representable by a linear differential equation 


with constant coefficients. The block diagram of a system 
having an input x and a corresponding output y is shown in 
figure 1. The symbol F(p) inside the box represents the 
system operator which, acting on the input x, yields the 
output y. The linear differential equation relating x and y 
may be written as 

d n y . d n ~ 1 y . . 

an dt n + a ’ 1 - 1 dt a - 1 ^ +a ° y 

, d m x . , d m ~ l x , . , 

— " m dt m ' ® m_1 dt m ~ l ‘ (1) 

If the system is initially at rest (or in equilibrium), term- 
wise application to equation (1) of the Laplace transformation 
(reference 2, pp. 51-56), defined as 


W)]=J o “/(«)e- p ‘ it 


and factoring of the result give 

(a n p n +a n .. 1 p’ t - l + . . . + a 0 )L(y ) 

=(b m p m +b m - l p m ~ 1 +. . .4 -b a )L(x) 
or 


T M b ™P m +b m -iP m ~ 1 +--- + b a 
yy a, n p n +a n -iP n ~ 1 -{- . . . +a 0 

Equation (3) may be formally written as 

L{y)=F(p)L(x) 


L(x) 


where 


F(ri) b ™P m + b m-iP m ~ 1 +- • -+6 0 
^ a„ p n +a n - l p a -'-+ . . . +a 0 


( 2 ) 


( 3 ) 

( 4 ) 


The transfer function of the system F(p) is defined as the 
ratio of the Laplace transform of any normal response of the 
system to the Laplace transform of the input producing that 
response. Normal response is the response of the system 
when initially’' at rest (reference 2, p. 26). 


X 


P(P) 


*- 


y 


Figure l.— Block diagram of linear system. 


From equations (2) and (4), the system transfer function 
may be written as 

m L{x) 


or 


| ye~& dt 

F(p)^ 

xe~P‘ dt 

EQUILIBRIUM FINAL CONDITIONS 


( 5 ) 


Derivation of frequency-response function. — The frequency- 
response function F(ia) is formally obtained directly from 
equation (5) by replacing p by iu (reference 1, pp. 96-98). 
Then, for a sinusoidal input x of frequency w, the output y 
ultimately is sinusoidal at the same frequency but with a 
relative magnitude and phase angle equal to the magni- 
tude and the phase angle of the complex number F(ia>). 
The frequency-response function is defined as 


F(ia) 



ye~ i( *‘ dt 


xe~ iu ‘ dt 


( 6 ) 


The term e~ iat is oscillatory and the integrals of equation (6) 

may not converge unless |y| dt and \x\dt converge. 

These conditions unduly restrict the choice of functions 
that may be used in equation (6) because of the implied 
requirement that x and y vanish as t increases without limit. 




FREQUENCY RESPONSE OE LINEAR SYSTEMS FROM TRANSIENT DATA 


549 


As will be shown, the previous restriction on x and y may 
be removed by suitable modification of equation (6). A 
much wider selection of input and output functions is thus 
permitted. 

For a system initially at rest, 

L[f'(t)]=pL[fm 00 


An alternate expression for the transfer function then is 


F(p)= 


W) 

L(x') 


( 8 ) 


or, from equations (2) and (8), 



y'e~ vt dt 


x'e pt dt 


that followed in -this report in connection with equation (10). 
The method of reference 4, however, differs from that of this 
report in being restricted to functions that begin and end 
at zero. In addition, the concepts of a system transfer or a 
system frequenc-y-response function are untreated. 

Modifications of frequency-response derivation.' — The 
frequency response may, at times, be desired directly in 
terms of the input and output functions themselves rather 
than their time derivatives. For example, if an electronic, 
device were used for analysis, the avoidance of the added 
complexity and attendant inaccuracy of differentiating cir- 
cuits might be desirable. 

If both x and y come to a steady state by or before some 
time t f , equation (9) may be rewritten as 

f f y'e~ u dt 

F(iu)= 4^7 — (ID 

I x'e~ ,a ‘ t dt 


The frequency-response function then becomes 

| y'e~ iut dt 
F(iu)^ Jo 




- ,u dt 


( 9 ) 


when p is replaced by ia. 

For any input or output function ending in equilibrium, 
x'~*Q and y'-*0, respectively, for The integrals of 

equation (9) therefore converge. 

From Euler’s relation, 

e -t“t_ cos o it — i sin at 


because for t > t r the integrands of equation (.11) are zero. 
Integrating equation (11) by parts gives 

rt f . . t f 

ice ye ~ lU dt-\-ye ,ui 
F(i«)- Jo -1 


or 


icjJ* f xe~ iu dt-\-xe w | 

(J^cos t ddt—y f t *V sin at dt+y f 0 ^^ 

^cos at dt — x rf' ) — i ^ J" 'ssin utdt+Xf 00 ^^ 

( 12 ) 




equation (9) may be written 


f* CO Z 1 CQ 

y' cos utdt—i y' sin at dt 

TV ; N */0 V 0 

X ^ rZ — 

x' cos ut dt—i x’ sin ut dt 
Jo Jo 


( 10 ) 


This form of the frequency-response function can be used 
with any data for which the input and the output begin and 
end in a steady state. 

The integrals of equation (10) can be evaluated in various 
ways. One theoretically exact method, which has been used 
at the NACA Lewis laboratory, utilizes a rolling-sphere 
harmonic analyzer. This type of analyzer, the same in 
operating principle as the device described in reference 3, 
was designed to obtain the Fourier coefficients from cyclic 
data, but the manner in which the coefficients are determined 
results in evaluation of integrals of exactly the form of those 
in equation (10) . The frequency-response function can thus 
be obtained for any frequency by operating the analyzer over 
the output curve to evaluate the integrals in the numerator 
and over the input curve to determine the integrals in the 
denominator of equation (10) without replotting the data in 
any other form. 


The frequency spectra for nonperiodic phenomena are de- 
termined in reference 4 by evaluating Fourier integrals with 
a rolling-sphere harmonic analyzer in a manner similar to 


An alternate procedure, leading to a somewhat simpler 
expression, employs the principle of superposition. For 
example, y may be considered the linear combination of a 
step function whose ordinate is y r and a second function 




(a) System-output time function. 

(b) Components of system-outpnt time function. 

Fiuuee 2.— Separation of system-output time function into components for equilibrium final 

conditions. 
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whose ordinates are those of y displaced downward by the 
value of y f . Such a separation of y into its component 
functions is shown in figure 2. An analytical expression of 
the process may be derived as follows: The transfer function 
of the system is 

I ye~ pt dt 

(5) 

J xe~ pi dt 

Because of the linearity of the transformation 

f ( y—yM~ p ‘dt+ f y,e~ pt dt 

— - — % 

(x—x f )e~ pt dt+ j x f e~ p ‘dt 


f (y-Vf)e~ pt dt+^ 
.Jo p 

f\x-x,)e-»'dt + ! 


Setting p—io} then gives 




(a) Piecewise continuous time function. 

(b) Piecewise approximation to time function. 

Figure 3.— General approximation to time function. 


Equations (12) and (13) define F(iu) for all frequencies 
except w=0 and F( 0) can readily be shown to be yjlzf. 

Approximation to frequency-response function. — Let /(£) 
be a piecewise continuous function as shown in figure 3. 
If within each of the r+1 pieces /(t) is approximated by a 
polynomial of any degree not greater than g— 1, then, for 
. t r+ i at infinity, the Laplace transform of J(t) is given by 

L\jm=± i Sk (f*+ o) — (t h - o)] (A io) 

k — o j= i y 

Equation (A10) is developed in detail in appendix A. 

The expression for L(f) in equation (AlO) requires only 
the values of the approximation and its derivatives at f=+0 
[/(+0), /'(+0), and so forth] together with the jumps 
(4+0) — / (s_Ji) (f* — 0)] in/ and in the derivatives of the 
approximating polynomials at any points of discontinuity 
(4) of the functions or their derivatives. The degree of the 
approximating polynomials may he different in the various 
pieces but g— 1 is the highest degree of polynomial in any 
piece. 

As an example of the use of equation (AlO), let/(£) be a 
continuous function (except at f=0) approximated, as shown 
in figure 4, by a series of straight lines. Here the poly- 

Tima function 


Figure 4.— Linear approximation to time function. 

nomials are of first degree and g=2. Hence, equation (AlO) 
becomes 

r 2 p—Vh 

/(/)= g g ^ L/' 2 -* (4+ 0) -/«-* (4- 0)] (14) 

Because / is continuous except at £=0 

/(4+ 0) -/&- 0) = 0 for A> 0 

/ (4 + 0) — /(4— 0) =/ (+ 0) for A— 0 
and 

iO)=p{ge-^[/ , fe+0)-/(< s -0)]+/'(+0)}+i/(+0) 

For evenly spaced t’s 

4+i — 4=4 

/'(4+ 0)=^±P^ 
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and for 0 


with 


Then 


/fc~Q)= /fc /* -1 

H 

Uf'(tr J rO)=f T -i-l — /r=0 

l / i ~{ (+Q - ) -+p/(+Q)] 


=4-1 S e-^</ t+ i+/ l -i-2/d+L/i-/(+ 0)Rpfi/(+ 0) 

p^l ( fc=l 


(15) 


or 


F{p) 


JM 0 i W) 

' £(x0 "*■ i(*0 

L(yi)+pL(y 2 ) 
£(x') 



(18) 


In order to get an expression that converges for p=0, 
multiply equation (15) by p and use the relation 


L(f)=pL(f) 


(7) 


Then, on setting p=ica and separating real and imaginary 
parts as in equation (9), equations (16) and (17) are obtained, 
in which the integrals converge for co=0. 


j f cos of dt=^ [2/*— (it-i+/* + i)] sin cakt i+/(+ 0) (1 6) 
I* f sin cat df=”| L/i-/(+0)]— S L+Zi+OJcos <akt x j 


where 


kti — f* 


(17) 


By applying equations (16) and (17) to the input and output 
curves of any system and using the results in equation (10), 
the approximate frequency-response function of that system 
may be determined. 

OSCILLATORY FINAL CONDITIONS 

The procedures thus far described apply only when the 
system comes to rest at some new equilibrium condition. 
If the output does not reach equilibrium but continues to 
oscillate about some mean value, a modification in the 
method is necessary in order to account for the oscillatory 
component. 

System linearity permits resolution of the output y into 
two components yi and y 2 , whose algebraic sum is the original 
output, as shown in figure 5. The block diagram equivalent 
to this resolution of y into components is given in figure 6. 
The dashed outline indicates the over-all transfer function of 
figure 1. From figure 6, F(p) may be represented as the 
linear sum of two transfer functions, F x (p) and F 2 (p), operat- 
ing in parallel. The following equation may then be written: 

F(p)=F 1 (p)+F 2 (p) 



(a) System-output time function. 

(b) Components of system-output time function. 

Figure 5.— Separation of systemoutput time fnnction into components for oscillatory final 

conditions. 


I 

l 

l 

X ^ 

Fi(P) 

+ Vl 

t 

t 



1 

1 

^ r 
1 
1 
1 



[ 

1 


F,(P) 

V* 

1 



1 

1 




'•--Ftp) 


Figure 6.— Block-diagram representation of linear system for oscillatory final conditions. 

If the steady oscillation is simple harmonic, the oscillatory 
component of the output may be expressed as 


y 2 =A sin (/3f-f <f>) 


for f>0. Then 


pL(y 2 )= 


Ap 


and 


P 2 +i S ; 


(/3 cos 4>-}-p sin <£>) 


Ap 


\ y x e~ v ‘ dt+ JY R - 2 (0 cos 6-bp sin <j>) 
F(p)=— (19) 


/; 


x'e~ pt dt 


When p is replaced by ica in equation (19), the frequency- 
response function becomes 


F(iu)= 


y x ' cos edit — sin Vi sin cat dt — cos 


( x' cos catdt—i ) x' sin catdt 
Jo Jo 


( 20 ) 
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TBEATMENT OF DEAD TIME 


EXAMPLES 


Dead time in a linear system is the time difference be- 
tween the initiation of a disturbance to the system and the 
beginning of a response to that disturbance. The mathe- 
matical equivalent of dead time then is a translation of the 
output function along the time axis in a positive direction; 
that is, if 

L[yffl=F(p)L[x® ] (4) 

for a system without dead time, then for the system with 
dead time 

L[y(t-M)]=G(p)Z[x(t)] (21) 


where G(j?) is different from F(p). The relation between 
G(p) and F(p) is readily determined by an expansion in a 
Taylor's series of y(t—At) in the neighborhood of t: 


y(t-At)=y®-My'(t)+^ v"'®+ ... 

Then 


For initial conditions of equilibrium 


L[y»m=P n L[y®\ 

Hence, 


=e~ Atp L[y(fy] (22) 

From equations (21) and (22) 

L[y(t)]=e^G(p)L[x(t)] (23) 

Comparison of equations (4) and (23) shows that 

e Atp G(p)=F(p) 

or 

G(p>)—e~ Alp F(p) 


from which the corresponding frequency-response functional 
relation is 

G(ia) = e~ M “F(io>) (24) 

It can easily be shown that the effect of dead time on the 
frequency-response function is to decrease the phase angle 
algebraically by the amount A tea, leaving the amplitude ratio 
unchanged. In general, F(iu) is a complex number and may 
be expressed as 

F(i(>))=Re iS 

From equation (24) 


G(ioi)=e iAf “Re ie 
or 

G(io})=Re^-^ 


Examples of the application of the method of this report 
to determine the frequency-response function of several 
systems are found in figures 7 to 9. The use of the- function 
thus found has not been covered because the detail with 
which such use has been treated elsewhere (for example, in 
references 1 and 2) makes extensive treatment in this 
report unwarranted. 

Examples 1 and 2 illustrate the types of frequency response 
obtained from several different types of time function. 
A method of this report was used with the specified time 
functions to determine the points shown in figures 7 and 8. 
Frequency-response functions determined directly from the 
known transfer functions of the systems used were also 
plotted for comparison. Example 3 presents data for a 
turbine-propeller engine and shows the type of frequency 
response obtained, using a method of this report. 

EXAMPLE 1 


A second-order system having the transfer function 


F{ ^ p 2 +2fp+l 


(25) 


is shown in figure 7 for both the underdamped and critically 
damped cases. The frequency-response function of the 
system is obtained directly from equation (25) by replacing 

pbyioi. For f 


F(iu)= 


1 

— w a +ici>+ 1 


o> 


and for f =1 


F(iu)= 


1 — o) tJ rbF 1 — c^+w' 1 

1 


(26) 


— 2iu~ f - 1 

1 — 0! 2 . 2 CO 

-% ■ 


(1 + co 2 ) 2 ■'(1+co 2 ) 2 


(27) 


System response to a unit step input is given for f =~ by 

A 


and for f = l by 


» = 1 _ ^ sin (~r <+ f) 


y=l— 


(28) 

(29) 


Frequency-response curves calculated from equations (26) 
and (27) are shown in figures 7 (c) and 7 (d). The rolling- 
sphere analyzer was used on the transient-response curves 
determined from equations (28) and (29) to obtain from 
equation (10) the points shown on the frequency-response 
curves of figure 7. Because a unit step input was used, the 
denominator of the right-hand side of equation (10) is equal 
to 1 for all values of frequency and F(io>) can be obtained 
by operation of the analyzer over the output curves only. 


EXAMPLE 2 


A first-order system is shown in figure 8. If the system 
has no dead time, the transfer function is 
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TT5 

and the corresponding frequency-response function is 

1 . T03 


1+tV 1+rV 
For a unit step input the output is 


(31) 


l=l-e T 


(32) 


If the system has a dead time At, the transfer function be- 
comes 

p-Atp 

(33) 

and 

G(io})=e- lAtu F(iu ) (34) 

The response to a unit step input then is 

t- At 

y=l-e T (35) 

for At. The frequency-response curves in figures 8 (d) 
and 8 (e) were calculated from equations (31) and (34). 
The function F(ica) was found from equation (31) and G(iu) 
was then determined from F(iu>) by increasing the phase 
angle of F(iu) negatively by the amount Atu, leaving |.F(i«)| 
unchanged. The rolling-sphere analyzer was used on the 
output curves of equations (32) and (35) to obtain from 
equation (10) the points shown on the frequency-response 
curves of figure 8. The use of a unit step input to the system 
makes the denominator of the right-hand side of equation 
(10) equal to 1 for all values of frequency. The frequency- 
response function can therefore be obtained by operation of 
the analyzer over the output curves only. 

EXAMPLE 3 

Experimental transient data obtained at this laboratory 
for the speed response of a turbine-propeller engine to changes 
in propeller-blade angle are shown in figure 9 (a). Because 
the frequency-response function is expressed as a ratio, 
units of the variables involved are unimportant and the 
ordinate of the data curves has been calibrated in relative 
values, the difference between initial and final values of the 
curves being set equal to unity. 

The frequency-response function F(i<S) of this engine for 
response of engine speed to changes in propeller-blade angle 
was found from equation (10). The rolling-sphere analyzer 
was used on the data. In figure 9 (b), F(iu) is shown as a 
frequency-locus plot in the complex plane; figure 9 (c) is a 
plot of amplitude ratio and phase angle against frequency. 






(a) Input and output data for turbine-propeller engine. 

(b) Frequency-response function of turbine-propeller engine plotted in complex pltmo. 
(e) Amplituderatioandpbaseangleofturblne-propellerfrequeney-responsefunetionagainst 

Figure 9.— Frequency-response function of turbine-propeller engino. 
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The closeness with which the data approach a semicircle in 
the complex plane indicates that the response of the system 
is substantially first-order. Tins inference is further borne 
out by the amplitude-ratio plot. The asymptotes and the 
solid curve were calculated for an exact first-order system. 

SUMMARY OF RESULTS 

Methods are presented for the determination of the 
frequency-response function of any linear system from 
general correlative time-response data for input and output 
of that system. Equations were developed for system 
equilibrium and oscillatory final conditions and for the effect 
of system dead time on the frequency-response function. 


For illustrative purposes, a method of this report was 
applied to several linear systems for which the frequency- 
response functions were known. 

Experimental data obtained at this laboratory for the 
speed response of a turbine-propeller engine to change in 
propeller-blade angle are presented and analyzed. 


Lewis Flight Propulsion Laboratory, 

National Advisory Committee for Aeronautics, 
Cleveland, Ohio, April 1, 194-9. 


APPENDIX A 

DEVELOPMENT OF EXPRESSION FOR GENERAL APPROXIMATION TO FREQUENCY-RESPONSE FUNCTION 

The Laplace transform of / may be defined as 

2*/)- [y*~ pt dt 


Let /(f) be a piecewise continuous time function approxi- 
mated as shown in figure 3. Thus, 


where 


/(f) = S M) 

/*© =u(t — 4-0«(4~ f) 9&) 


(Al) 


(A 2) 


and u(t— t k ~\)u(t k — t) is a pulse of unit height and duration 
fit— 4-i- Use of the pulse for defining in equation (A 2) 
insures vanishing of the function outside the interval 
4 -l<£< 4- The polynomial g k (t) may be of any degree 
not greater than g— 1 and may be of different degree in 
different intervals. 

As shown in appendix B, the gth derivative of /*(£) is 
A- (5) (f)=u(t— 4-iM4— t)ff* <e) (0+2 (f— 4- 1 )£4 C4-;f) (4-0 — 

>=i 


x u™(t-t k )g k «-»(Q 
l-i 

Because g k (t) is of degree <g— 1, 

<7/ 5l (£)=0 


(A3) 


Then from 


and equation (A3) 


/«>(i)=g /*«© 


Jt=*l ’ 


(A4) 


(AS) 


f a) (0=2 2 li® (i-4-i)sr*"-' ) (4-0 - 

fc=i j=i 

/ Ca) ®=EE "(jt-tt (4) - (4)3 + 

fc=l j=l 


2 v/»(t-t r+1 )g T+l «-n(t r +i) (A6) 

J=1 7=1 


(A7) 


Then, under the assumed zero conditions at t— — 0, 


Lu w (M=p«Lm\ 

L [u (i> (t — 4)] =p 7 ’~ l e~ rC t 


and, as f r+l -»-oo, 
Let 




(A8) 


L[u U) (t~t r+ d\-*0 

(4) (4+ 0) 


ff t+ i' e - 1) (4)=/ le - fl (4+0)) 

> (A 9) 

0fc C *~ fl (4)=/ c *- fl (4-O) ) 


with/ Ci /5 (— 0)=0 imderstood. Then equation (A6), when 
combined with equations (A8) and (A9), becomes 

i[/(0]=2 2 [/<*"* (4+0) -/*-* (4- 0)1 (A10) 

i=0 7=1 P 

which is the equation used in the text. This expression for 
L[f(f)] includes the values of /(£) and its derivatives at 
£= + 0 but not at f=— 0. 

Note that if the Laplace transform of / had been defined 
as 


£(/)= P/<r p ‘rff 

J+o 


(All) 


the same answer would have been obtained for L(j) with 
— 0)=0 understood whether such is the case or not. 
This result necessarily follows because / is not impulsive at 
f=0. When equation (All) is used, 


LU (1) m=P Q L[f(t)\ - ± (-L- 0) (Al 2) 


and 


7=1 


L\u if) (f)I = 0 for j> 1 


(A13) 
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APPENDIX B 

DEVELOPMENT OF EXPRESSION FOR /*'«>«) 


In appendix A is defined by 

jk(t) =u(t—t k -i)u(t k —t)g k {t) (A2) 

Assume, for the moment, that the derivatives of a finitely 
discontinuous function such as the step and an infinitely 
discontinuous function such as the impulse have meaning. 
Then, a formal differentiation of equation (A2) gives 

u'(t— 4-iM 4— t)0s©— 4-i}w'(4— (B 1) 

The impulse u' is different from zero only at its discontinuity. 
Hence, 

4-i)«(4— t)g K '{t)+u'(t— 4-i)ff*(4-i)— ' u'(t k — t)g k (t k ) 

(B2) 

because 

m( 4 — 4-i)=l 

Because u' is an even function 

u'(t k ~t)=u'(t—Q 
and equation (B2) becomes 

ft® =u(t—t k . i)«(4 ~ t) g k ' (t) + u’ (f — 4_j) gr*(4- 1 ) —u'(t— 4) <7*(4) 

(B3) 

Differentiation of equation (B3) gives 

=u(t — 4_i)u(4— t) g/ (i) +u'{t—t k ^)u{t k —t)g k '(t) — 

u(t — 4- iK (4 — (t) + u"(t — 4- i)ff*(4-i) — u'% — 4)0t(4) 

(B4) 

If the reasoning used on equations (Bl) and (B2) is applied 
to equation (B4), 

ft"® =u(l— 4-iM4— t)g k "(t)+ [u’(t~ 4-i)0*'(4-i) + 

v!'($— 4-i)fl , s (4-i)] — [«' (f — 4)?/(4) 4)^(4)] 

CB5) 

Continuation of this process shows finally that, in general, 
(f) =u(t-4- 1 )«(4-<)^* (9) © +g « (fl 0-4-i)^ Cs -«(4-i)- 

(A3) 

1=1 

An heuristic justification of the procedure for obtaining 
equation (A3) will now be given. 

A commonly used method of defining the single impulse is 
to begin with a pulse such as that used in equation (A2) and 
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then, allow the pulse duration to approach zero while main- 
taining constant the enclosed area of the pulse (reference 5) . 
Doublet, triplet, and higher-order impulses may be defined 
similarly but such a means of definition does not lend itself 
to a consistent interpretation of the meaning of the deriva- 
tives of an impulse. 

A better means is found in the use of an infinitely differ- 
entiable function to define the pulse. Then, the derivatives 
of the function are the derivatives of the pulse and when, in 
the limit, the pulse becomes the impulse, the derivatives of 
the pulse become the derivatives of the impulse. A conven- 
ient function is the normalized Gaussian error distribution, 
defined by 



h(t—t k )= — \=.e 2 

V. »• / 

(B6) 


vy2x 



From 





1 

\ 

1 

II 

1 — 1 

for all <r 

(B7) 

and 

t/ — 




lim h(t— f*i)=0 




< r -»0 



lim k(t— Q = -r 00 
0--> 0 




k{t—t k ) is seen to satisfy the requirements for an impulse 
occurring at t~t k . Equation (B7) holds for a— »0; hence 
h(t—tk ) satisfies the requirement that the step be its integral 
as or — >0 or, conversely, that 

lim 

c-id 


By successive differentiations, followed by letting ff-»0, it 
can be seen that 


f*)=lim h'(t—t £ 

o —> 0 

u"%— 4)=lim h"{t—t k ) 

4 T — >0 

and, in general, 


f*)=]im t k ) 

o — *0 


The function h{t— t h ) and its first three derivatives are shown 
in figure 10. Note that h and its even derivatives are even 
functions; the odd derivatives are odd functions. 
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